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Abstract 

We describe the hydrodynamic behavior of the fc-step exclusion pro- 
cess. Since the flux appearing in the hydrodynamic equation for this 
particle system is neither convex nor concave, the set of possible solutions 
include in addition to entropic shocks and continuous solutions those with 
contact discontinuities. We finish with a limit theorem for the tagged par- 
ticle. 



1 Introduction and Notation 

In his paper Liggett (1980) introduced a Feller non conservative approximation 
of the long range exclusion process to study the latter. A conservative version of 
this dynamics, called fc-step exclusion process was defined and studied in Guiol 
(1999). It is described in the following way. 

Let k e N* := {1, 2, ...}, X :— {0, 1}^ be the state space, and let {XnjneN 
be a Markov chain on Z with transition matrix p(., .) and P^{Xo — x) — 1. 
Under the mild hypothesis sup^^g^ v) ^ Lk, defined below, is an 

infinitesimal pregenerator: For all cylinder function /, 

Lkfiv) = E <ik{x,y..v)[f{v''''')~fiv)], (1) 

where qk{x, y, rj) = JTiZi ^ ^(^Oj 5: , CTj, < is the intensity for mov- 
ing from a; to y on configuration r/, (jy — inf {n > I : Xn — y} is the first (non 
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zero) arrival time to site y of the chain starting at site x and r/^'" is configuration 
T] where the states of sites x and y were exchanged. 

In words if a particle at site x wants to jump it may go to the first empty 
site cncoimtcrcd before returning to site x following the chain Xn (starting at 
x) provided it takes less than k attempts; otherwise the movement is cancelled. 

By Hille-Yosida's theorem, the closure of generates a continuous Markov 
semi-group Sk{t) on C(X) which corresponds to the fc-step exclusion process 
{Vt)t>o- Notice that when k = 1, {f]t)t>o is the simple exclusion process. An 
important property of fc-step exclusion is that it is an attractive process. 

Let Xfe be the set of invariant measures for {rit)t>Q and let S be the set of 
translation invariant measures on X. lip[x, y) = p{0, y — x) for all x, y G Z and 

.) is irreducible then 

(Jfen5)e = {!^a:ae[0,l]}, 

where the index e mean extremal and Ua is the Bernoulli product measure with 
constant density a, i.e. the measure with marginal 

Va{ri e X : r]{x) = 1} = a. 

In this paper we prove conservation of local equilibrium for the totally asym- 
metric process in the Riemann case i.e.: p{x,x + 1) = 1 for all x & Z and the 
initial distribution is a product measure with densities A to the left of the origin 
and p to its right, we denote it by ii\,p. The derived equation involves a flux 
which is neither concave nor convex and appears for the first time as a hydro- 
dynamic limit of an interacting particle system. Up to now the "constructive" 
proofs for hydrodynamics relied on the concavity of the flux, see Andjcl & Varcs 
(1987) or the papers by Seppalainen {e.g. Seppalainen (1998)), whose key tool 
is the Lax-Hopf formula. The entropy solution for the type of equation we con- 
sider was first studied in Ballou (1970). Such solutions can have entropy shocks 
as well as contact discontinuities. Our aim is to take advantage of Ballou's 
result to deduce conservation of local equilibrium also in a constructive way, 
in the spirit of Andjcl & Varcs (1987). However we explain in the last section 
how to derive hydrodynamics for general initial profiles, and nearest neighbor 
dynamics. 

In section 3 we prove a law of large numbers for a tagged particle in a fc-step 
exclusion process. 

2 The hydrodynamic equation 
2.1 Heuristic derivation of the equation 

Since the process is (totally) asymmetric we take Eulcr scaling. For every r G R 
define ?7f (r) := rj^-ifde^^r])^ where [e^-'^r] is the integer part of e^^r, for e > 0. 

Given a continuous function u^{x),x G R (initial density profile), we define 
a family of Bernoulli product measures {t'^o}£>o on X, by: For all x G Z, 

i^'^o{veX:7]{x) = l} = u°{ex). 
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We call {j^^o} the family of measures determined by the profile u". Let u^{r, t) := 
J{Sk{£~^t)r]Q{r))di'^o{'rio)- Then for all r S M, M^(r, 0) converges to u°{r) when 
e goes to 0; applying the generator (|l|) to r]^{r) we have 



-5fe(te-i)(^g(r)) = £-i5,(te-i) 



fc-i 



-J2I['^oir + je) [1 -Vo(r+i^ + 

i=0 j=Q 



k-1 i 

Assuming that local equilibrium is preserved (thus expectation of products fac- 
tor), and taking expectations with respect to v'^o^ we obtain 



fe-l i 
■i=0 i=0 



If we now let e converge to 0, u(r, t) lim^^o u^(r, t) should satisfy 



du dGk{u) 



dt 



dx 







(2) 



u(x, 0) — 

where represents the flux of particles: 

k 

Gfc(u) = ^jV(l-u). 

This is a non standard form because Gk is neither convex nor concave, thus 
equation (^) is no longer "a genuinely nonlinear conservation law" , using the 
language of Lax (1973). To deal with this equation, we have to use an extended 
version of non linear Cauchy problems treated by Ballou (1970). 

Remark 2.1 

Let k go to infinity and denote by Goo the limiting flux function: 

Goo(m) = -r~~- 

That case corresponds to the totally asymmetric long range exclusion process. 
The resulting equation is simpler because the flux function Goo is strictly convex. 
The hydrodynamics in this case should follow from the arguments of Aldous & 
Diaconis (1995) for the Hammersley's process. 

For notational simplicity, from now on we restrict ourselves to the case k = 2. 
However our arguments can be easily extended for all k. 
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2.2 Hydrodynamics in the Riemann case 
2.2.1 Notation and result 

Our main theorem characterizes the hydrodynamic (Euler) limit of the 2-step 
exclusion process at points of continuity, when the family of initial measures is 
determined by a step function profile. From the heuristic derivation we would 
expect that in the hydrodynamic limit the density profile would satisfy equation 
(H). We show that the limiting density profile at time t is the entropy solution 
of equation (||) starting with the initial value a step function profile. We 
now give a brief summary of results concerning the solution of equation 
due to D.P. Ballou (1970), when u° is a step function. This will motivate the 
formulation of the theorem as well as some aspects of the proof. 

Existence of weak solution to the Cauchy problem given by equation (||) 
with bounded measurable initial condition was proved in Ballou (1970), under 
the assumptions: 

1. GkeC^R). 

2. G'j^ vanishes at a finite number of points and changes sign at these points. 
In order to obtain uniqueness further conditions are needed. We require our 

solutions to satisfy the: 
Condition E: (O.A. Olemik) 

Let x{t) be any curve of discontinuity of the weak solution u{t,x), and let 
V be any number lying between :— u{t,x{t) — 0) and m+ := u{t,x{t) + 0). 
Then except possibly for a Rnite number of t, 



S[v;u^] > S[ 



+ . 



^ _ u 
where 

Gk{w)-Gk{v) 

b\v]w\ :— . 

w — V 

It is known (Ballou (1970)) that the following two conditions are necessary 
and sufficient for a piecewise smooth function u{x^ t) to be a weak solution of 
equation (||): 

1. u{x,t) solves equation (|^) at points of smoothness. 

2.1ix{t) is a curve of discontinuity of the solution then the Rankine-Hugoniot 
condition (i.e. d[x{t))/dt = 5[u+;u^]) holds along x{t). 

Moreover condition E is sufficient to ensure the uniqueness of piecewise smooth 
solutions, which are the entropy solutions to the equation. Hereafter we only 
deal with the case k = 2. We denote G{u) := G2{u). 

If G were convex (concave) only two types of solutions would be possible. 
We now describe these two types of solutions. 
Let u"{x) = Al{^<o} + pl{x>o}- 

If A > p (p > A), then the speed of characteristics which start from a; < 
(given by G") is greater than speed of characteristics which start from a; > 0. If 
the intersection of characteristics occurs along a curve x{t), then since 

+ _ GjX) - G{p) 
S[u^;u ] = ^ = S[\] p] 
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Rankine-Hugoniot condition will be satisfied if x'{t) = S[X; p]. Thus 

"^"^'^^-t p, x>S[X;p]t. 

is a weak solution. The convexity of G implies that condition E is satisfied 
across x{t). Therefore u{x,t) defined above is the unique entropic solution in 
this case and will be referred to as a shock in the sense of Lax (1973). 
If A < p (p < A), then the characteristics starting respectively from .t < 
and from x > never meet. Moreover they never enter the space-time wedge 
between lines x = Xt and x = pt. We can choose values in this region to obtain 
a continuous solution, the so-called continuous solution with a rarefaction fan: 
Let h be the inverse of H := G' , 

( X, x< H{X)t; 

u{x,t) = l h{x/t), H{X)t < X < H{p)t; 
[ p, H{p)t<x. 

It is possible to define piecewise smooth weak solutions with a jump occurring 

in the wedge satisfying the Rankine-Hugoniot condition. But the convexity of 
G prevents such solutions to satisfy condition E. Thus the continuous solution 
with a rarefaction fan is the unique entropic solution in this case. 

For the 2-step exclusion process, the flux function G is neither concave nor 
convex. Instead G{u) = u + v? — 2u^ is convex for m < 1/6 and concave for 
u > 1/6. In this case in addition to the shock and continuous solution with 
a rarefaction fan it is possible to have solutions for which the ciirve of discon- 
tinuities never enters the region of intersecting characteristics. The quotation 
in boldface is the original number of lemmas, prop... in Ballou (1970), but the 
notation refers to 2-step exclusion: 

Definition 2.2 [B def2.1] 

For any u < 1/6, define u* := u*{u) as 

u* = sup{r] > u : S[u; r]] > S[v; u] Vu e {u, r])}. 
For any u> 1/6, define := «*(«) as 

u* = inf{?7 < u : S[u-,ri] > S[v;u] Vf G {r],u)}. 

In other words, for u < 1/6, if we consider the upper convex envelope G"^ 
of G on (u, +oo), then u* is the first point where G° coincides with G. In the 
same way when u > 1/6, is the first point where the lower convex envelope 
Gc of G on {—oo,u) coincides with G. For rj < 1/6, 77* = (1 — 277)/4, and for 
V > 1/6, r?. = (1 - 27?)/4. 

Let hi and /12 be the inverses of H respectively restricted to (—00, 1/6) and 
to (1/6, -hoo). I.e. hi{x) = (1/6)(1- V7- 6a;) and h^ix) = (1/6)(1 -h V7 - 6a;) 
for X G (—00, 7/6). 

The following lemmas are taken from Ballou (1970). 
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Lemma 2.3 [B lem2.2] Letrj < 1/6 be given, and suppose thati]* < oo. Then 
S[v;v*]=H{r,*). 



Lemma 2.4 [B lem2.4] Let t] < 1/6 be given, and suppose that if < oo. Then 
rj* is the only zero of S[u\ rj\ — H{u), u > rj. 

If A < p < 1/6, the relevant part of the flux function is convex and the 
unique entropic weak solution is the continuous solution with a rarefaction fan. 

If p < A < p* (p < 1/6), then H{r]) > H{p) if p < 77 < 1/6, and H{r]) > 
H{p*) > H{p) if 1/6 < rj < p* since H is decreasing in this region. Thus 
H{X) > H{p), which implies an intersection of characteristics: The unique 
entropic weak solution is the shock. 



Let p < p* < X {p < 1/6): Lemma 2.4 applied to p suggests that a jump from 



p* to p along the line x = H{p*)t will satisfy the Rankine-Hugoniot condition. 
Since p* is specially defined for this, a solution with such a jump will also satisfy 
condition E. Therefore if we can construct a solution with the jump described 
above it will be the unique entropic weak solution in this case. Notice that 
since H{X) < H{p*), no characteristics intersect along the hue of discontinuity 
a; = H{p*)t. We call this case contact discontinuity, following Ballon. The 
solution is defined by 

\, x< H{X)t; 

u{x, t) = { h2{x/t), H{X)t <x< H{p*)t; 
9, H{p*)t < X. 

Corresponding cases on the concave side of G are treated similarly. 

Let T denote the shift operator. We are able now to state our result. 

Theorem 2.5 Let u G M, X, p ^ 1/6, and pA,p the Bernoulli product measure 
on TL with densities X for a; < and p for a: > 0. Then 

lim p\.pTi^t]S2it) = iyu(v,i) 

at every continuity point o/m(.,1), where iyu(v,i) denotes the product measure 
with density u{v, 1) defined by: 

Case 1. X < p < 1/6; continuous solution, with a rarefaction fan 

r A, x< H{X): 

u{x,l)=< hi{x), H{X) < X < H{p); 
[ p, H{p) < X. 

Case 2. p < X < p* , (p < 1/6): entropy shock 
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Case 3. p < p* < X, (p < l/Sj: contact discontinuity 

( A, x<H{\)- 
u{x,l) = l h2{x), H{\) <x<H{p*)- 
{ p, H{p*) < X. 

Case 4- 1/6 < p < A; continuous solution, with a rarefaction fan 

( A, x<HiX); 
uix,l) = l h2{x), H{X) <x<H{p); 
[ p, H{p) < x. 

Case 5. p > X> p», (p > 1/6): entropy shock 

^ I P, x>S[X;p]. 

Case 6. p > p^, > X, (p > 1/6 j; contact discontinuity 

A, x<H{X); 
u{x,l) = { hi{x), H{X) <x<H{p,); 
p, H{p^) < X. 



Remark 2.6 

For any k > 2 the profiles will be of the same kind, because Gk has only one 
inflection point between and 1 and is first convex then concave. 

Remark 2.7 

Comparing with hydrodynamics of simple exclusion we observe that fc-step ex- 
clusion (fc > 2) has not only a stable increasing shock (Case 5) and a decreasing 
continuous solution (Case 1) but also a stable decreasing shock (Case 2), an 
increasing continuous solution (Case 4) and two contact discontinuities (Cases 
3 and 6). 



2.2.2 Proof of Theorem 2.5 



It follows the scheme introduced in Andjel & Vares (1987), where the authors 
obtained the hydrodynamic limit for the one-dimensional zero-range process in 
the Riemann case, i.e. the hydrodynamic equation 

du dchiu) 

h = 

dt dx 

U{X, 0) = U°{x) = Al{a;<o} + pl{x>0} 

was derived. There (j), the mean flux of particles through the origin, was a 
concave function. Therefore, their proof used both the monotonicity of the 
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process, still valid here, and the concavity of the flux, that we have to replace 
by an ad hoc use of the properties of the solution of (H) . 

Informally speaking, they first showed that a weak Cesaro limit of (the mea- 
sure of) the process is an invariant and translation invariant measure. Then they 
showed that the (Cesaro) limiting density inside a macroscopic box is equal to 
the difference of the edge values of a flux function. These propositions were 
based on monotonicity, and on the characterization of invariant and translation 
invariant measures (both valid for k-step as well) ,thus we can quote them (with 
appropriate notation for the fc-step), and take them for granted. 

Lemma 2.8 [AV 3.1] Let ji be a probability measure on {0, 1}^ such that 

(a) Vp < fi < i^x for some < p < X < I, (b) either /iri < or fiTi > fi. 
Then any sequence Tn oo has a subsequence Tn^. for which there exists D 
dense (countable) subset ofM. such that for each v G D, 



for some /it, G X2 H 5. 

Lemma 2.9 [AV 3.2] For v (z D, we can write fx^ — J Vajvida), where is 
a probability on [p, A]. Also, if u < v are in D, 



with, for w E D, F{w) = j[wa — G(a)]'^y,{da). 

The difficult part is then to prove that 7„ is in fact the Dirac measure 
concentrated on ^(a;, 1). They did it in Lemma [AV 3.3] and Theorem [AV 2.10] 
using the concavity of their flux function. 

For fc-step exclusion, at this point we will have to look separately at the six 
cases given in the theorem. 

To conclude, Andjel & Vares had to prove that the Cesaro limit implies a 
weak limit, through the following propositions, based on monotonicity. We will 
use these results also without proof. 

Proposition 2.10 [AV 3.4] Let /i = fix,p- If 





) 



F{v) - F{u) 



(3) 




ifveD,v< S[X;p] 
ifveD,v> S[X;p] 



then 
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Proposition 2.11 [AV 3.5] If n satisfies 

(a) ^ < vx, (h) fiTi > fi, (c) there exists vq finite so that 

1 

lim - / ^iT[^t]S2{t)dt = vx 

for all V > Vq. Then 

lim fiT!^t] S2 [t) — Vx for all v > vq. 



Proof of Theorem 2.5 in Cases 2 and 3. In case 1, the proof is not different 
from the one given in Andjel & Vares (1987) since G is convex in the relevant 
region; thus we omit it. In case 2, G is not convex in the relevant region, 
therefore we supply a proof, though it is quite close to the original one. Case 
3 uses ideas from cases 1 and 2 and introduces some new ideas to deal with 
complications arising from the non-convexity of G. Cases 4-6 are symmetric 
to cases 1-3 in the sense that the roles played by convexity and concavity are 
exchanged. 

The proof has 3 steps. The main ingredients are monotonicity, and inequal- 
ities relying on the properties of S[.; .], G, H given before. 
First step. 

Using monotonicity of the 2-step exclusion, we can proceed as in the begin- 
ning of the proof of Lemma [AV3.3], and get two finite values v and v so that: 
If w e D and v > v, then = Sp, while — Sx H v < y. 

Second step, preliminary. 

Let u < V, both in D. Attractiveness of the process (since Vp < fix,p < ^x) 
and (^) imply 



{v — u)p<v J a^y{da) — J G{a)^y{da) — u J a^uida) + J G{a)^u{da) 

<{v- u)\ 

(4) 

(i) Taking u < y, the first step gives — so the second inequality of (||) 
is simplified in 

V J ajy{da) — J G{a)^y{da) ~ uX + G{\) < {v — u)\ 

which can be written 



(G(A) - G{a))jy{da) <v (A - a)-fy{da) (5) 
(ii) Similarly, for v < v, ^y = 5p, and the first inequality of (|4|) reads 

(a - p)7„(da) < / (G(a)-G(p))7„(da) (6) 

[P,A] J[p,\\ 
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Proof for Case 2: p < A < p*, p < 1/6. 

The definition of p* implies that for every a € (p, A) 

G{\) - G{a) ^ GjX) ~ G{p) 
A — a A — p 

so that inequahty (ph for v < w < S\\\ p] {v G D) yields 



(A - a)7„(da) > S'[A;p] / (A - q;)7i,(q;). 
b.A] hp A 

Since v < S[X; p] we conclude that = Sx- 

Starting from inequality (j^) with v > u > S[X; p] (u e D) and proceeding in 
a similar manner we can show that "fu — 6p. 

Proof for case 3: p < p* < A, p < 1/6. 
Second step, part 1. 

Let u < V, both in D, u < y, and v < H{X}. On [p*, A), G is concave, thus 
for every a E [p* , A] 

H{X) <S[a:X]<H{p*). (7) 
If a e [p, p*) then from the definition of p* it follows that 

S[a;p*]>H{p*). (8) 

We decompose [p, A] = [p, p*] U (p*, A], and (^) becomes 

(A — p*)7t,(da) + w / (p* — a)7„(da) + u / (A — a)7„((ia) > 

(G(A)-G(a))7,(da)+ / (G(A) - G(a))7,(da). 
[p,p*] -^(p'.A] 

By i) 

(G(A) ~ G(a))7„(da) - / [(G(A) - G(p*) + G(p*) - G(a)] 7„(da) 

[P,P*] "'[P:P*] 

> / [(A-p*)5[p*;A]+i/(p*)(p*-a)]7.(da) 
-'[P^P*] 

and by (|) 

/ (G(A) - G(a))7„(dQ) > F(A) / (A - a)-/y{da). 
Those two inequalities together with (||) give 

{v - S[p*- A])(A - p*) f ^y{da) + {v~ H{p*)) f (p* - a)7„(da) 

"'[P:P*] "'[P,P*] 
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+ {v-H{X)) {X~a)-/vida)>0. 
Since v < H{X), by the only possibility is 7„ = (5a- 

Second step, part 2. 

Let u,v G D, such that v > v, and H{p*) < u < v. We decompose each 
integral of (||) on the two intervals [p, p*] and (p*, A]. 



{G{a) ^ G{p)hu{da) 



(P*,A] 



(G(a) - G(p*))7„(da) + / iG{p*) - G{p)hu{da) 

(P*,A] "'(P%A] 



By definition of p* and by Lemma p.3| , for all a G [/?, p*], 
G(a)-G(p)<i/(p*)(a-p), 
and G being strictly concave on {p* , A] , we have for all a G [p* , A] 
G(a)-G(p*)<i7(p*)(a-p*). 

Thus by (|) 

u/ {a - p)ju{da) + u {a - p*)^u{da) + u j [p* - p)^u{da) 

J[p,p'] J(P':>-] J{P':>-] 

<[ H{p*){a- p)-iu{da)+ I H{p*)ia - p*)-f^{da) 
Ap.p'] "'(p%a] 



which can be written 



iu-Hip*)) 



+ H{p*){p* ~ p)^u{da) 

■J(p*M 



[a - p)^u{da) + I [a - p*)^u{da) 
[p,p-] J(p'^] 



< 0. 



+ / (P* - P)lu{da) 

'(P*,A] 



Since u — H{p*) > 0, this implies 



[P.A] 



(a - p)^u{da) = 



Therefore we conclude 7„ = 5p. 
Second step, conclusion. 
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Using the two preceding parts, attractivity, Propositions 2.10 and 2.11, we 
conclude in case 3 

and in case 2 

lim M.,,r[„,5,W = ( (11) 
t-»oo [ J \ ' l_ z^p, 1/ u > A[A; ^ ' 

Third step, first part. 

Let ui, w, wi e -D, ui < i?(A) < d < H{p*) < vi, such that ui and wi belong 
to an interval where G is concave. By attractivity, 

fl \ 

limsup/XA,p52(t) - V < (« - "i)A (12) 

V N*l / 

The second step and (||) imply 

/ [v^t] \ 

\\m^ix^pS2it) - r^{x) \^vip- G{p) - mX + G(A) 

\ [uit] J 

which, combined with ( p^ ) gives 

V [vt] J 

Since ^^(A) < v < H{p*), there exists some 9, p* < 9 < \, with v — H{9), 
so that 9 — h2{v). Let 9 < 9' < \, we now apply ( p^ to /ig'.p, and we use 
attractivity through /ig'.p < /iA.p to get 



(13) 



/ ^ [fit] \ A bit] 

hminf ^A,pS'2(t) - V ?7(a^) > hminf ^(,,,p52(0 - V ?7(a;) 

\ ["f] / V / (14) 



>vip~G{p) + G{9')-v9' 

By (0) again. 



AKt] \ 
limmf MA,pS'2(t) -^5I'?(^) - ^^iP-G(p)-'y / 
V [vt] J -^I' 

which, together with (|l^), il wc make ^ 0, gives 



ajv{da)+ / G(Q;)7^((ia) 
[p,A] J[p,\] 



G{9) -v9< (G(a) - va)'y^{da) (15) 

-'[p,A] 
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Since G{a) is convex when a G (p, 1/6) and concave when a £ (1/6, A), for 
all w, G{a) — va has at most two critical points determined by the condition 
G'{a) — V. One satisfies 1/6 < a < A and is a local maximum, the other (when 
it exists) is a local minimum with p < a < 1/6. We want to conclude that the 
local maximum is a global maximum when p < a < X. We know from Lemma 
{^) that 

Gip*) ~ G{p) = {p* - p)H{p*) 
which implies, since v = H{9) < H{p*), 

G{p*) ~ VP* > G{p) - vp 

Because G(a) — va is increasing in {p*,9) (recall that H{9) ~ v) we have: For 

aU a e {p*,e) 

G{a) -va> G{p*) - vp* > G{p) - vp 
We conclude that is a global maximum: 

max [G(a) - va] = G{0) - v0 

p<a<X 

thus = (5e = Sh2(v)- 

Third step, second part. 

This part follows closely the argument in Andjel & Vares (1987), but we 
detail it for the sake of completeness. Since the measures ^ /g P'\,pT[yt]S2{t) dt 
depend monotonically on v and form a relatively compact set, we have for all v, 

T 

t^\,pT[vt\S2{t) dt = Vh2(v)- 

It remains to prove that 

^}^t^\,pT[vt]S2{t) = Vh2(v) 

when H{X) < v < H{p*) (by continuity of /12, this result will also be valid at 
H{X)). For this, let p.y be a weak limit of fJ.x^pT[yt]S2{t)- It is enough to show 

(a)/i„ > I'h^iv), {b)fiv{v{0)) ^h2{v). 

(a) Let 9 = hiiv), with H{X) <v < H{p*), v = H{9). Let p* < 9 < 9 < X, 
then V < H(9) and p§ pT[vt]S2{t) < lJ-x,p'T[vt]S2{t)- Added to (|o|), this yields 

lim p.Q ^vi\S2{t) ^Vg< lim px.pTy^i\S2{t) = i^v 
Hence, by continuity, if 9 converges to 9, 

fly > V0 — l^h2 (v) ■ 



lim — / 

T^oo T in 
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(b) Let Ul < H{X) < H{p*) < vi. By the definition of u{x,t) in that case, 
u{x, 1) dx — X dx + h2{x) dx + p dx 



m Jh(\) JH{p') 

= X{H{X) ~ Ul) + J' 9H'{0) dO + p{vi ~ Hip*)) 

= vip - uiX - G{p) + G{X) 

where we have integrated by parts the integral of the second line (derived by a 
change of variables); therefore, by (||), 

(i \ n 

lim px.pS2{t) - ^ ri{x) = / u{x, 1) dx (16) 



and (b) follows from ( |16D and (a), with the same argument as in g] p. 332 (proof 
of Theorem 3.2), based on monotonicity. ■ 



3 Asymptotic behavior of a tagged particle 

We introduce here an interpretation of the /c-step exclusion dynamics valid in 
the totally asymmetric case. Up to now we considered that a particle might 
jump from x to the first empty site in {a: + 1, a; + k}. If we want to leave the 
particles ordered we could equally say that the particle at x pushes the "pack" 
of (< fc) neighboring particles in front of it, each one moving of one unit to the 
right. In other words, if there is no particle at site a; + 1 then the particle at 
site X goes to site x + 1 ; if there is one particle at site x + 1 and no particle at 
site X + 2 then the particle at x pushes the particle of site a; + 1 to site x + 2 
and occupies site a; + 1; and so on... Then the generator reads 

fc-l i 

LkRri) = E E n '^(^ + (1 - '?(^ + * + [/(??^^^+i'-^+'+i) - /(t?)] 

a;eZ i=0 i=0 

(17) 

where 

{■q{xi) ifu = a;i, 
7,(x,_i) iiu = x,, i = 2,...,l, 
rjiu) otherwise. 

It is easy to see, comparing ( |l7|) to (|l]), that they do correspond in this 
setting, because we do not label the particles. The interest of ( pT|) w.r.t. (|j) is 
that it keeps track of the particles' order. We will need this interpretation in 
the next section. 

Similarly, we define a Tagged "Pushing" Particle, and the generator of the 
fc-step exclusion process as seen from this Tagged Pushing Particle is 
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Lkiv) 



k 71—1 

+ E n vim){l-vin))[f{nv°-'-n-fiv)] 

n—l ni—l 

k—1 k—7i —1 m — 1 

+ E E n n - ^(")) - fiv)] ■ 



71—1 l — l m—~l i—1 

To be clearer, let us write and comment it for k — 2. 

L2iv)= E rK.x)(l-r;(:r + l))[/(r;---+i) -/(,?)] (18) 

0:5^0,-1 

+ (l-r7(l))[/(ri,yO-i)-/(,7)] (19) 

+ V{^)v{^ + m-v{^ + 2))[f{v^^'^'^^^')-f{v)] (20) 

x^0,-l,-2 

+ V{l){l-rm[f{riv"^''')~fiv)] (21) 

+ - ry(l)) [/(riry-i^°-i) - /(ry)] . (22) 

Part ( p^ ) is simple exclusion involving sites away from the origin, part (19) 
corresponds to the "classical" tagged particle for simple exclusion. Part (20) 
is a "strictly" 2 steps exclusion involving sites away from the origin. Part (21) 
describes the "pushing" of the tagged particle. Finally in part ( p^ ) the tagged 
particle is "pushed" by another particle. 

A straightforward adaptation of the simple exclusion case (see Ferrari (1986)) 
gives 

Theorem 3.1 The Palm measure Va of (i.e. the measure on X defined by 
Vai-) — '^ai-lvi^) — ^)) invariant and ergodic for the k-step exclusion process 
as seen from a tagged pushed particle. 

Sketch of proof: For instance when /c = 2 it is enough to show 

^^M-l)vml - vii))fiv-'^°-')} = i^MOMm - vm)f{Tiv°^'-'')} 

and 

M{?7(-2)r7(-l)(l - vmfiv-'^-''')} = M{rK-l)^(0)(l - v{l))f{riV-'''^')} 

which are obvious for any G 5 (recall that S is the set of translation invariant 
measures on X). ■ 
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Theorem 3.2 Law of large numbers for the Tagged Pushing Particle (k — 2). 

For a 2-step exclusion process with initial distribution v^, if Y{t) denotes 
the position at time t of a Tagged Pushing Particle starting at the origin then 

Y(t) 

lim — ^ = (1 - a)(l + 2a) Pp^a.s. 

t^oc t 

Proof: Using the notation of Ferrari (1992a) pp. 41-43 we define the instanta- 
neous increment of the position of the Tagged Pushing Partiele by 

m - lim nnt + h) Yit)\Y,^.) 
= (1 - r]{x + 1)) + r]{x - 1)(1 - r]{x + 1)) + r]{x + 1)(1 - t]{x + 2)). 

So 

■0 dVa = 1 — a + 2q;(1 — a) 

Furthermore limh^oE^^{Y {t + h) — Y{t))'^/h < +oo so that the conditions of 
theorem 9.2 of Ferrari (1992a) are satisfied, which gives the result. ■ 

Remark 3.3 

Referring to results of Guiol (1999), we notice that a Tagged Pushing Particle Xt 
behaves as a regular tagged particle in the fc-step exclusion. Indeed, intuitively, 
the "regular" tagged particle can make long jumps, so is expected to move faster, 
but it cannot be pushed; and the rate at which the Tagged Pushing Particle 
moves compensates exactly those long jumps. 



4 Generalizations 

1. We obtained conservation of local equilibrium for the Riemann case. It can be 
generalized to an initial product measure with a non-constant profile of bounded 
variation following the papers of Rezakhanlou (1991) and Landim (1993) (see 
also Kipnis & Landim (1999) chapters 8 and 9). Indeed Condition E is equivalent 
to Kruzkov entropy inequality (see Godlewski & Raviart(1991), Lemma 6.1 p. 88) 
which is the key tool of the latter proofs. Moreover, since we are working with 
a one-dimensional totally asymmetric process, we can equivalently consider its 
interpretation defined by (|l^); this way, particles do not jump over each other. 
We can therefore follow Section 6 of Rezakhanlou (1991) which contains a two- 
block estimate valid only for one-dimensional nearest neighbor processes, and 
enables to conclude the derivation of hydrodynamic equation without using 
Young measures. The key point consists in coupling two versions of the process, 
to prove that the number of sign changes between them can only decrease, 
due to attractivity and the reduction to a one-dimensional nearest neighbor 
case. For fc-step exclusion, this result is a straightforward adaptation of Liggett 
(1976), Lemma 5.1. For the same reason, extension to deterministic initial 
configurations is possible, using Venkatsubramani (1995). 
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2. We treated only totally asymmetric /c-stcp exclusion but our proof is also 
valid for nearest neighbor asymmetric transition rates, which lead to a flux 
function with the same properties. However one has to compute carefully this 

function, due to the complexity of the fc-stcp dynamics. For instance for fc = 5 
and transition rates p{x, x + 1) = p, p{x, x — 1) = q = 1 — p, p > q, we obtain 

G|'«(u) = «,(! - u) [{p -q){{l + 2u) + 3u^il-pq) 

+4v^{l - 2pq) + bu'^il - 3pq + p^q^)} + 3u'^2p^q] . 

Indeed the last term corresponds to a "cycle" jump: To go from a; to a; + 3 when 
sites X + 1 and x + 2 are occupied, the particle follows the path {x,x + l,x + 
2,x+l,x + 2,x + 3). 

The extension to an initial product measure with a general profile is still 
possible, using Rezakhanlou (1991) and Landim (1993), but with the help of 
Young measures. 
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